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Abstract: SHiP is a newly proposed fixed-target experiment at the CERN SPS with the aim 
of searching for hidden particles that interact very weakly with Standard Model particles. The 
work presented in this document investigates SHiP’s physics reach in the parameter space of the 
Neutrino Minimal Standard Model (vMSM), a theory that could solve most problems unexplained 
by the Standard Model by incorporating sterile neutrinos. A model introducing an extra t/( 1) 
symmetry in the hidden sector, providing a natural candidate for dark matter, is also explored. This 
work shows that the SHiP experiment can improve the sensitivity to Heavy Neutral Leptons below 
2 GeV by several orders of magnitude, scanning a large part of the parameter space below the B 
meson mass. The remainder of the vMSM parameter space, dominated by right-handed neutrinos 
with masses above 2 GeV, can be explored at a future collider. Similarly, SHiP can greatly 
improve present constraints on 17(1) dark photons. 
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1. The SHiP experiment and its physics case 

SHiP is a newly proposed general purpose fixed target facility at the CERN SPS accelerator. A 
400 GeV proton beam will be dumped on a heavy target in order to produce 2 x 10^° proton-target 
interactions in 5 years. A dedicated detector downstream of the target will be used to probe a 
variety of models with light long-lived exotic particles with masses below ^(10 GeV/c^). Active 
neutrino cross-sections and angular distributions will also be studied, thanks to a dedicated detector 
placed between the target and the hidden sector detector QE 0, ED- 

SHiP’s flagship goal is to use decays of charm and beauty mesons to search for Heavy Neutral 
Leptons (HNLs), which are right-handed partners of the Standard Model (SM) neutrinos. The ex¬ 
istence of such particles is strongly motivated by theory, as they can simultaneously solve multiple 
questions unanswered by the SM. In the Neutrino Minimal Standard Model (vMSM), HNLs allow 
to explain the baryon asymmetry of the Universe, account for the pattern of neutrino masses and 
oscillations and provide a dark matter candidate [0]. 

2. Estimating SHiP’s sensitivity 
2.1 Heavy Neutral Leptons 

Using the algorithms described in detail in |0], a Monte Carlo simulation was developed to deter¬ 
mine both the rate of HNLs produced at the target and the acceptance of the HNL decay products. 
Prom these estimates, the expected number of events in 5 years of SHiP operation was calculated. 

The number of expected HNLs can be parametrised by factorizing the problem into the HNL 
production rate and detector acceptance: 

n{HNL) = A(p.o.t) x x{pp ^ tlNL) x x s^tot{HNL visible) (2.1) 

where A(p.o.t) = 2 x 10^° is the number of proton-target collisions integrated over five years of 
operation; X{PP HNL) is the production rate of HNLs per pp interaction; is the probability 
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that an HNL of given mass and couplings decays in SHiP’s acceptance; and s/tot {HNL —visible) 
is the fraction of HNLs decaying in the detector acceptance and resulting in a detectable final state. 

The production cross sections of c- and b- mesons, the main sources of sterile neutrinos at 
SHiP, are estimated with PythiaS [01, and the HNL production rate takes into account all the 
dominant, kinematically-allowed decay channels of and mesons into sterile neurinos. The 
widths of these channels are parametrized according to [0 as a function of the HNL mass and 
couplings. The HNL lifetime is estimated as the sum of the widths of its dominant decay channels 
(3v, TT^V, 71^^, p°v, and final slates), lhal are also paramelrised according lo [0. 

PythiaS is used lo build a dala sel conlaining kinemalic informalion from chaimed and 
beauly mesons produced al Ihe SHiP largel. For every meson, all Ihe kinemalically allowed de¬ 
cays info massive sterile neulrinos are simulated using Ihe ROOT TGenPhaseSpace class [0]. 
A fwo-dimensional binned probabilily density function (PDF) is built out of the resulting HNLs 
momentum p and polar angle 6, weighted with the branching ratio of the meson decay in which 
the HNL is produced. This PDF was then used as an estimate of the 4-momentum spectrum of 
HNLs produced at the SHiP target. Each bin of this PDF contributes to the total vertex acceptance 
proportionally to its content and to the integral of the decay probability over SHiP’s decay 
volume, as a function of p and 6: 

^v.x(p,0)= (2.2) 

Jorilr 

^vtx = j S^vtx (p, 0) dpde (2.3) 

The detectable fraction of HNLs is computed as a combination of the branching ratios and 
final sfafe accepfances of fhe HNL decay channels fhaf are kinemalically allowed and produce fwo 
charged parficles in fhe final slate: 

j^tot{HNL visible) = ^BR(//AL ^ i) x £/{i) (2.4) 

i 

where fhe index i runs over visible final slates and each acceplance (/) is compuled by simulation. 

The HNL yield depends on fhe hierarchy of fhe active neulrino masses, fhe HNL couplings fo 
fhe Ihree SM flavours Ug,U^,U^ and fhe HNL mass. The fofal coupling Uf and fhe HNL 

mass are free parameters in fhe simulation, while fhe relalive sfrenghfs of fhe coupling and fhe hier¬ 
archy are fixed and consfilute differenl scenarios which conform lo existing Iheorefical sfudies [0. 
Figure |I] shows SHiP’s sensilivify, along wifh fhe foreseen sensifivify of FCC-ee in fwo realis¬ 
tic configuralions (Section in fhe paramefer space of fhe vMSM, for a scenario wifh inverfed 
hierarchy of acfive neulrino masses and relalive coupling slrenglhs 17^ : ~ 48 : 1 : 1. 

2.2 Dark Photons 

A variely of New Physics models foresee a sel of SM-neufral unobserved parficles fhaf do nof 
inferacl wifh SM parficles excepl fhrough a “messenger” parlicle belonging fo fhe hidden secfor. 
Minimalistic models consisf of a 17(1)' gauge symmefry in fhe hidden secfor, whose gauge boson 
/ is called a dark photon. If 17(1)' is broken by a Higgs-like mechanism, fhe dark photon can 
acquire a non-zero mass. Such dark photons may be searched for af SHiP in neulral di-lepfon and 
di-meson final slates: / may mix fo fhe SM pholon fhrough loops of parficles charged bolh under 
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Figure 1. Physics reach in the vMSM for SHiP and 
two realistic FCC-ee configurations (see text). Pre¬ 
vious searches are shown (dashed lines), as well as 
the cosmological boundaries of the model (greyed- 
out areas) [Ein- 



Figure 2. SHiP sensitivity to dark photons produced 
in proton bremmstrahlung and secondary mesons de¬ 
cays. Previous searches [HU explored the greyed- 
out area. Low-coupling regions are excluded by Big 
Bang Nucleosynthesis and SN1987A UTO, 0. 


(7(1) and (7(1)^ Assuming / is the lightest particle of the hidden sector, it would then decay to 
and qq final states thi'ough a virtual photon. 

A toy simulation was used to estimate SHiP’s sensitivity to dark photons. The algorithm is 
analogous to the one described in Section o and to the method used by the authors of [Hi. 

Figure @ shows SHiP’s sensitivity to dark photons, compared to previous searches. The con¬ 
straints from supernovae cooling [0 and big bang nucleosynthesis, and the most up-to-date limits 
established by previous experiments up to July 2014 IHDU . are shown in grey. 

3. HNL searches at future colliders 

A review of possible methods to perform HNL searches at future colliders is given in [imi . 
This document reports on the sensitivity that can be achieved in direct searches in Z decays. FCC- 
ee is a 90-400 GeV high-precision collider aimed at precision tests of the Standard Model, 
currently studied within the scope of the Future Circular Collider project at CERN. Luminosity 
studies suggest that the machine could be operated at the Z resonance to produce 10*^ to 10^^ Z 
bosons per year WTSl . 

A method similar to the one outlined in Section o was used to compute the expected number 
of recorded HNLs. HNL production is assumed to happen in Z —s- vv decays with one neutrino 
kinematically mixing to an HNL. If the accelerator operates at the Z resonance, Z bosons decay at 
rest and the HNL lifetime is boosted by a factor Y = mz/2m[^ + mj^/2mz. All final stales are 

considered detecfable. A detector with spherical symmetry is modeled. Backgrounds from W*W*, 
Z*Z* and Z*y* processes can be suppressed by requiring a displaced secondary vertex. 

The minimum and maximum displacements of the secondary vertex in FCC-ee, referred to as 
r in Figure |T], depend on the resolution and dimensions of the tracking system. The two FCC-ee 
sensitivity contours shown in Figure |I] represent two realistic scenarios HTTll . The production of 
lO'^ (lO'^) Z bosons is assumed. Z bosons are allowed to decay in v +HNL pairs, followed by the 
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subsequent HNL —)■ Uv decay. Secondary vertices between 1 cm (100 /im) and 1 m (5 m) from 
the interaction point are selected. 

This work shows that the SHiP experiment can improve the current limits on Heavy Neutral 
Leptons by several orders of magnitude, scanning a large part of the parameter space below the 
B meson mass. Similarly, SHiP can greatly improve present constraints on dark photons. Right- 
handed neutrinos with larger mass can be searched for at a future Z factory. The synergy between 
SHiP and a future Z factory would allow the exploration of most of the vMSM parameter space for 
sterile neutrinos. 
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